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We discuss the exact solution for the properties of the recently introduced "necklace" model for 
„Cj ■ reptation. The solution gives the drift velocity, diffusion constant and renewal time for asymptoti- 

cally long chains. Its properties are also related to a special case of the Rubinstein-Duke model in 
Mh , one dimension. 
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■ I. INTRODUCTION 

o 

The most common lattice model for reptating polymers is the Rubinstein-Duke model Q, Q . The major drawback 
of this model is, that it has sofar defied an exact solution. Therefore it is natural to look for models that are in 
the same spirit, but which are more open to analytic rather than numerical methods. Recently Guidoni et al. Q 
introduced and analyzed a 1-dimensional model for reptation, to which we shall refer as the "necklace model" . It is a 
chain that moves through the exchange of beads and vacancies along a line. The purpose of this note is to show that 
the model permits an exact solution by means of the Matrix Product Expansion, due to Derrida et al. 0. Also the 
spectrum of the Master Operator can be derived, leading an exact expression for the renewal time. In addition we 
relate the properties with those of a special 1-dimensional variant of the Rubinstein-Duke model with hernia creation 
, and annihilation. 

a 

o . 

C->, II. THE NECKLACE MODEL 



The necklace model is a string of N + 1 beads located on a line of points. The beads are either neighbors or nearest 
neighbors. In the latter case there is an unoccupied lattice site (a vacancy) between the beads. The beads are not 
allowed to occupy the same lattice point and two consecutive vacancies are forbidden in order to ensure the integrity 
of the chain. The possibility of vacancies is an expression of the elasticity of the chain. The internal beads hop with by 
, exchanging with a vacancy. The two end-beads can also exchange with a vacancy from outside. Clearly the number 
of beads is conserved, but not the number of vacancies, which enter and leave at the ends and migrate through the 
chain. The stationary state of the system follows from the Master Equation, governing this stochastic motion. 

The key to the exact solution is to focus on the motion of the vacancies rather than on that of the beads. There are 
N positions for the vacancies available, since each vacancy must be surrounded by beads. Each of the N positions can 
be occupied by a vacancy or not. The state of the chain of beads is fully determined by the occupation distribution 
of the vacancies. So we consider the chain as an open system for vacancies, which hop with the rates derived from the 
motion of the original beads. Their hopping is constrained by the rule that two vacancies can not occupy the same 
position. It would mean that the two beads surrounding such a double vacancy would be separated by two vacancies. 
Transfering the motion for the beads to the vacancies we arrive at the following rules 



• An internal vacancy can hop to a neighboring empty position with rate p c to the right and p' c to the left. 

• The exchange of a vacancy with outside goes with rates p a for entering and p' a for leaving at the left hand end 
of the chain. 

• On the right hand end, the rates are pb for leaving and p' b for entering. 

Note that we have primed rates for the hops to the left and unprimed for the motion to the right. The asymmetry 
between left and right can be attributed to a driving field on the beads, inducing a bias B for a hop to the right and 
a bias B~ l for a leftward hop. These biases may derive from a charge on the beads, which is influenced by an electric 
field. Then we have for the ratios the relations 

— — ^ — ^ — b 2 (1) 

Pa Pb Pc 
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For the solubility of the model it is not neccessary to assume these ratios and only when we discuss the properties in 
more detail, we will use this physical restriction. The above given rates are more general than those used by Guidoni 
et al., in particular we allow for a finite driving field. For our exact solution, it is however important that the hopping 
rules for the internal vacancies are uniform along the chain. 



III. THE MASTER EQUATION 



The state of the chain can be represented by N variables r<, which assume the values 1 and 0, where 1 viz. 
corresponds to the presence viz. absense of a vacancy. The probability distribution P(r l5 • ■ • , t^) of the stationary 
state follows from the Master Equation 

MP(n,---,T N )=0, (2) 

with M. the Master Operator. It contains the usual gain term consisting of all transitions which increase the probability 
P(ti, ■ ■ ■ , Tj\r) and a loss term which contains all transitions out of the configurations (n, • • • , tjv). The Master Operator 
is the sum of N + 1 operators representing the action of the beads 

N 

M = Y J M l . (3) 

As first example consider the bead at the left hand of the chain, leading to A^o- It acts on the T\ dependence of 
P(ti, • ■ • , tat), as it only influences the existence of a vacancy on the first position. It is given by the expression 

M P(n ■••) = (!- rx)\p' a P{l ■■■)- Pa P(0 ■■■)}+ n\p a P(0 ■■■)- p' a P(l ■ ■ ■)}. (4) 

The first 2 terms refer to the case where a vacancy is absent from the first possible position. It has a gain and a loss 
term. Similarly the last 2 terms refer to the case where a vacancy is present. The expression can be shortened to 

MoPin ■ ■ •) = (l - 2 n ) [ P ' a P{\ ■ ■ ■) - Pa P(0 ■ ■ ■)]. (5) 

In the same way the other parts of the Master Operator are expressed as 

[ M : P(- ■ ■ T ]lT]+1 ■■■) = ( Tj - r j+1 )\p' c P(- • • 0, 1 • • •) - PcP(- • • 1, • • •)] 

i (6) 
[ MnP{---t n ) = (l-2T N )\p b P(---l)-p' b P(- ■■())}. 

This form of the Master Operator shows that it has an eigenvalue 0, as summation over all Ti yields zero for all 
individual terms, implying conservation of probability. The left eigenvector is constant for all configurations, but the 
right eigenvector is the non-trivial probability distribution of the stationary state. 



IV. MATRIX PRODUCT EXPANSION FOR THE STATIONARY STATE 



The idea of the matrix product expansion is to represent the probability distribution P(ti, ■ • ■ , tat) as product of 
matrices 

N 

P( Tl , ■■■,t n ) = Z~ 1 (W\ UfaV + (1 - t 3 )£]\V). (7) 

i=i 

As we shall see in a minute, T> and £ are sort of creation and annihilation operators, represented by infinite dimensional 
matrices and \V) and (W\ are states in this space. The normalization is given by 

Z N = (W\(V + £) N \V). (8) 

The matrices T> and £ , as well as the states (W\, and \V) are to be determined such that Q) is the stationary state 
of the Master Equation. By insertion of (J7J we find the action of the Mj on the factors 

MolnV + (l - n )£] = (i-2 n )(p' a V-p a £), 
Mj [tjV + (1 - tj)£] [r j+1 V + (1 - r j+l )£] = ( Tj - T J+1 ){p' c £V - Pc V£), (9) 
M n [t n V+{\~t n )£] = (l-2T N )(p b V-p' b £). 
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Then we use the identity 

(r 3 - r j+l ){V + £) = {!- 2T j+1 )[ Tj V + (1 - Tj )£\ - (1 - 2t 3 )[t 3+1 V + (1 - r j+1 )£] (10) 
and impose the condition 

p' c £V - p c V£ = £[D + £], (11) 

in order to rewrite the expression for the internal reptons as 

MjfcV + (1 - t 3 )£][t 3+1 V + (1 - Tj+1 )£] = 
C(l - 2r 3+1 ) [t 3 V + (1 - t 3 )£] - C(l - 2r,-) [r j+1 V + (1 - r J+1 )£]. 



(12) 



The "commutation" relation (jl 1|> defines the properties of the V and £ matrices. The states \V) and (W\ are fixed 
by the relations 

(W\(p a V-p a £) = (W\(, (p' b £-p b V)\V)=£\V) (13) 
Then it is possible to recombinc in (J5J) the factors again to probabilities 

( M P(t 1 ,---,t n ) = ((l-2n)P(T 2 ,---,T N ) 

M 3 P(n, - ■ ■ ,T N ) = -((l-2Tj)P(Tx,---,Tj-i,T j+ l,---,T N ) 

(14) 

+ C(l - 2t 3 + 1 )P(t 1 T 3 ,T j+ 2 7 - ■ ■ , T N ) 

K M N P(n,---,T N ) = -C(1-2t n )P(t 1 ,---,t N -i) 

The parameter £ is arbitrary; it has been introduced to facilitate the normalization. The sum over all these relations 
vanishes, which shows that (J7J) is the stationary state of the Master Equation. 
Relation (|ll(l can be clarified by the substitution 

V=l + b\ £ = l + b (15) 

This yields for the creation and annihilation operators W and b the g-deformed commutation relation 

btf - qtfb =l-q, (16) 

provided that we fix £ and q as 

C=p'c~Pc q=Pc/p' c - (17) 

The advantage of the new operators is that their spectrum has been thoroughly investigated. They have similar 
properties as the usual creation and annihilation operators for the harmonic oscillator. Note that for q = 1 the 
operators seem to commute, but it is easy to rescale w and b such that the right hand side of Hlfijl becomes 1, as often 
is done @- 

The power of the representation is that it gives the probability distribution for all chain lengths N . The properties 
of the matrices T> and £ and the states (W\ and | V) are independent of N. It requires, however, quite a bit of formal 
manipulations to retrieve the chain properties from the general expression (|7|. 

The relations (| 1 1 ft and (|13J) were mentioned by Derrida et al. jj. They were related to the g-deformed algebra by 
Sasamoto and Blythe et al. @ , leading to explicit formulae. Unfortunately, the case that we see as most physical 
(relations J3Jl and all intrinsic mobilites the same), does not fall in the wide range of parameters already treated. 
However, the technique that they employ, can be used to deal with our case. As it is not our aim to marginally extend 
the exact solutions to an even wider regime, we refer to 0,0 for details. 

V. THE DRIFT VELOCITY 

The phase diagram as function of the parameters is governed by the magnitudes of the net hopping rates 

Apa = Pa ~ Pa, &Pb = Pb ~ Pb, &Pc = ^(p'b ~ Pc) (18) 
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The smallest is the limiting factor. If the input on the left hand side is small, a dilute phase results. If the output at 
the right hand side is small, a dense phase will form. If the throughput in the bulk of the chain is the limiting factor, 
a maximum current phase appears, with a constant density of vacancies in the bulk. 

For illustration we discuss the case (0, with hopping rates only influenced by the driving field with bias B. Then 
we are in the maximum current phase. We restrict ourselves to the most interesting property: the drift velocity. The 
drift of bead j, in a specific configuration, is given by 

vj p' c P(-- -0,1 •••) /»,/';••• 1 . i) .. (19) 

The first term is the probability to jump to the right and the second for that to the left. Remember that the beads 
jump in the opposite direction of the vacancies. Comparing this with relation ©, we have the identity 

TjMjPin,- ■ ■ ,T N ) = Tj(l - T j+1 ) Vj. (20) 

By summing over the variables Tj and Ty+i, the factor in front of Vj gives 1 and one finds for the average drift 

( Vj ) = ( Tj Mj). (21) 
This average can be expressed, with the aid of l|14|) as 

v = Z N 1 (W\(V + S) N ~ 1 V) = ^(. (22) 

We have dropped the index j as the result is obviously independent of j. It has to, since there is no accumulation of 
beads in the stationary state. So l|22|l is the expression for the drift velocity, relating it to the normalization factors 
Zjq. The proportionality of the drift to £ = p' c — p c is simply a matter of times scales: when the rates go up, the drift 
follows proportionally. 

The Zn are expressed in terms of what is called the position operator + b of the q-deformed harmonic oscillator. 

V + £ = 2 + b + b j . (23) 
The spectrum of W + b has been investigated in detail . It can be represented as 

(b + b f )\9) = 2 cos 0\0) (24) 
The spectrum is continous: < 9 < ir. So we get 

Z N = [ d6w(6)[2(l + cos 6)} N (W\6)(9\V) (25) 
Jo 

Here w(9) is a function that enters in the closure relation for the states \9) (see |j|). For large N the power of H-cos6* 
gets strongly peaked around 9 = and one can apply the saddle point method to find the leading term. We put 9 = 
in all non-singular terms yielding 

Z' N ~ z(q)4 N N-V 2 (26) 

where the factor z(q) is independent of N and contains amongst others the factor (V7|0)(0|y). Fortunately this factor 
drops out in the drift velocity. 

v ~i(p' c -p c ) (27) 

The drift vanishes of course when the bias B — » 1, but its leading term does not decay with the power TV -1 as does 
the drift in the RD-model. 

The diffusion coefficient is usually defined in the limit that the driving field vanishes B — > 1. Using the Einstein 
relation one has to divide by a factor N and this leads to the (dimensionless) value 

D ~ — (28) 
AN V ' 

Relations l|27|) and l|28|l show that we are in the maximum current regime, which is not limited by the input and 
output at the ends of the chain. 
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VI. SMALL CHAINS 



The asymptotically leading behavior is of course the most interesting result, showing the strength of the Matrix 
Product Expansion. However, it is also worth while to see that the whole behavior is determined by the commutation 
relation (|16[1 . To illustrate this point we confine ourselves to the physical restriction (JTJ. Then the equations for (W\ 
and \V) reduce to 

(W\tf = q(W\b, b\V) = qtf\V). (29) 

Note that (W\ = (V\ since the first relation l|29(l is the conjugate of the second. 
As first example of playing with commutation relations consider 

{V\(b + tf)\V) =q{V\(bt+b)\V) (30) 

where we used ()29(l to interchange b and by applying b to \V) and b' to As q ^ 1 the average has to vanish. 
This holds also for all odd powers of b + b< . The second example is the average of b'b, for which we first use the 
relations 

(V\btf\V) =q- 2 (V\tfb\V) (31) 

But the commutation relation (|16J) can also be used to convert the right hand side of Ij31|l to the left hand side 
expression. From these two relations one deduces 

(V\tfb\V) = ^_~JI (V\V) (32) 

All other averages of 2 creation or annihilation operators directly follow, with (|29|l . from this expression. So the 
current of the 2 link system equals 

_ 2(1 + q + q 2 ) , 

VN = 2 - TT6q-TW { } (33) 

an expression, which of course, can also be obtained by solving the probability distribution from the Master Equation. 
By this technique it will be a long and hard road to get to the behavior at large N. 



VII. THE RENEWAL TIME 



Another interesting quantity is the renewal time. It is defined as the slowest time of decay towards the stationary 
state and it follows from the spectrum of the Master Operator as the eigenvalue with the smallest negative real part. 
Clearly the corresponding state decays the slowest and the renewal time is the inverse of the gap in the spectrum. 
The matrix product representation does not lead to the full spectrum of the Master Operator; it only gives the 
stationary state eigenfunction. On the other hand the renewal is usually defined in the fieldless case for which the 
Master operator becomes hermitian (or symmetric in our case) . The terms of the Master Operator can be expressed 
in terms of the operators Oj, creating a vacancy, and aj annihilating a vacancy (see Q). 

Mo = B(a\ - a\ai) + B^ 1 (a\ — aio{), 

Mj = B(a]a J+ i - Ojo]a] +1 Oj + i) + B^ 1 (a J a ] J+1 - a]ajOj+iat +1 ), (34) 
M = B~ 1 (a N - ajyaw) + B(a f N - ajvajy). 
The aj are hard core boson operators, equivalent with spin 1/2 operators, via the relations 

a f = (a x +ia v )/2, a = (a x - ia y )/2. (35) 

with the consequence 



a ] a = (1 +ct 2 )/2, aa f = (1 -cr z )/2, 



(36) 
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where the cr'a are the Pauli Matrices. The Master operator, expressed in these spin operators, is an unusual hamil- 
tonian since it is non-hermitian. However, for the undriven system B = 1, it turns into the well known ferromagnetic 
Heisenberg chain. 



N-l 

^=2El¥i+i-l]+< + 4-2. (37) 
j'=i 

Due to the magnetic field in the x direction on the boundaries of the chain, it is profitable to work in a basis of 
eigenstates of <j x . These are the symmetric (spin up) and the anti-symmetric (spin down) combinations of a vacancy 
and its absence. 

T> = 2- 1 / 2 [|1) + |0)], I i > = 2- 1 / 2 [ |1> - |0> ]. (38) 

The "groundstate" (being the highest in the spectrum) of the chain is the state with all spins directed in the x 
direction, which has an eigenvalue 0. In vacancy language this is the state in which all configurations have the same 
probability. 

As the hamiltonian l|37[) conserves the number of spins in the x direction the spectrum breaks up into sectors with 
a given number of spins up. The smallest excitation from the "groundstate" is in the sector with one spin down. Let 
x n be the value of the state with the down spin at the position n. Then we have the set of equations 

-(3 + A)^i + x 2 =0 

xi - {2 + \)x 2 + x 3 =0 

=0 (39) 
x N ^i — (3 + A)a;Ar = 

The eigenvalue spectrum is readily evaluated and one finds 

A(fc) = -2(1 -cos(7rfc/A0), k = 1, 2, • • • , N. (40) 

Thus the gap in the spectrum is given for k = 1 with 

A(l) = -2(1 - cos(7r/iV)). (41) 
The other branches can be evaluated similarly, but lead to larger negative eigenvalues. 



VIII. THE RUBINSTEIN-DUKE MODEL 



The Rubinstein-Duke (RD) model is in a way complementary to the necklace model. In the RD model the dynamical 
elements are reptons, blobs of monomers of the order of persistence length. A polymer chain is a string of reptons 
on a lattice with the constraint that reptons arc cither in the same cell or in neighboring cells. There are two types 
of links between reptons: the slack links connect two successive reptons in the same cell and the taut links those 
in neighboring cells. A repton can hop to a neigboring cell if it does not leave an empty cell behind. So the basic 
mechanism of motion is the interchange of taut and slack links. The elasticity of the chain is due to storing length 
through occupation of cells with more than one repton. The integrity of the chain is guaranteed by the requiring that 
a connected tube of cells is traced out by the chain. The natural embedding is in a lattice of dimension d > 1, but 
Duke 2] showed that it suffices to study the projection on the driving field direction. Thus the dimension d becomes 
a parameter and the links are characterised by a 3- valued variable: yi = for a slack link, yi = \ for a taut link in the 
direction of the field and yi = — 1 for a taut link opposite to the field. The parameter d influences the ratio of slack 
to taut links, but not the universal properties of the chain e.g. the exponents by which the renewal time and the drift 
velocity depend on the chain length. So the model is often studied for d = 1 although the 1-dimensional version is 
rather artificial 0. 

Sartoni and van Leeuwen || noticed that the 1-dimensional version of the RD model could be related to a model 
with two types of particles. Consider the superposition of two non-interacting 1-dimensional systems of particles. The 
particles are called + viz. — and they hop according to the same rules as the vacancies in the necklace model model. 
So they are not allowed to occupy the same site, but, since they are non-interacting, a + particle may occupy the 
same site as a — particle. Therefore each site can be in 4 states: empty, occupied by a + or a — or double occupied, 
which we indicate by a ±. In total there are 4^ configurations. As the two systems do not interact, the probabilities 
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of the combined system are the product of the probabilities of the two systems. Each of them can be treated with 
the Matrix Product Expansion and yields identical expressions for e.g. the drift velocity. 

Having the freedom to choose the transition rates of the two systems independently, we drive the + particles in one 
direction and the — particles in the opposite direction with equal strength. Then the + particles are identified with 
the links yi = 1 of the RD model and the — particles with the links yi = — 1. The and ± state are both mapped on 
a slack link = 0. This gives a contraction of the A N configurations of the superposition to the 3^ configurations of 
the RD model. Knowing the probability distribution of the combined system one can construct the transition rates 
in the Master Equation for the contracted system. 

It will in general be a quite involved calculation. We illustrate it by a simplifying assumption, which is motivated 
by the following consideration. By reversing the driving field, the role of particles and vacancies are interchanged. 
The system is invariant under a transformation that maps the particles of one component onto the vacancies of the 
other component. In this map an empty site, which is a vacancy state for both components, maps onto a doubly 
occupied site. Thus, as we will see, it is reasonable to equate the probabilities for an empty and a doubly occupied 
site. 

Before we list the result for the transition rates of this constructed RD model, we note that the underlying combined 
system has a feature not present in the usual RD models. If we have a succession of two taut links, it can be a 
combination of the and ± state. In that case, one of the particles of the ± state can move to the empty site, thereby 
creating a pair of opposite taut links. This is called the creation of a hernia. It is not incorporated in the usual RD 
model, since it is a modification of the tube (consisting of the set of taut links). Similarly the opposite process: a + 
and — particle meeting each other on the same site, is the annihilation of a hernia. With this in mind we come to the 
following list of transitions. 

• Transitions of a slack-taut combination. The slack link corresponds to the 2 states of the + and — particle 
system. Each of them can interchange the taut and slack link, so the transition rate is the same. 



Transitions of a slack-slack combination. This corresponds to 4 states of the particles. Two of them cannot 
move: the 0, and the ±, ± combination. The two others can create a hernia. So the hernia creation rate is 
1/2. 

Hernia annihilation. A hernia corresponds to a unique particle state and it may develop into a slack-slack pair 
in two ways: the + particle may move or the — particle. So the hernia annihilation rate is 2. 

A slack end-link. It has 2 particle configurations, and ±, and both may move to the taut position, either by 
creating a particle in the state or by annihilation of a particle in the state ±. So the transition rate is 1. 

A taut end-link. It is a unique particle state and it can transform itself in two ways in a slack state: by 
annihilation of the particle or by creation a particle of the other kind. So it has a transition rate 2. 



Thus the particle system maps onto a chain with asymmetric transition rates for the end reptons and for the creation 
and annihilation of hernias. The other moves are the same as in the RD model. The model with the above listed 
transition rules will demonstrate the same behavior as the superposition of the + and — particles system with opposite 
driving fields. Although the map is based on a global and not a local symmetry, it is correct in the weak driving field 
limit. We have calculated independently the diffusion coefficient and the gap of this model by means of the DMRG 
method 9]. Fig. ^ compares the drift velocity with the expression (|28f) . refined with finite size effects (as can be 
deduced from |8|]). 

N+l 

VN = (42) 

Note that the drift velocity of the RD model is twice the value of the necklace model, since there are two systems of 
particles moving independently. 

Fig. ^3 does the same for the renewal time and compares it with the expression (|4ip. As one observes the agreement 
is perfect. Looking into the numbers one has to conclude that the differences can be made arbitrary small by making 
the DMRG calculation more accurate. It proves also that the DMRG method is amazingly accurate in calculating 
reptating polymer chains. 



IX. DISCUSSION 



We have given an exact solution for the stationary state of the necklace model introduced by Guidoni et al. 0] • 
Also an exact expression is given for the low lying excitations (gap). This model mimics the reptation of a chain 
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FIG. 1: The drift velocity for the RD model with hernias as a function of the length of the chain, together with the expression 
1421 for the necklace model (a). Part (b) gives the renewal time of the RD chain with the expression (1411 deduced from the 
necklace model. 



in a 1-dimcnsional system. It is not difficult to formulate the model in higher dimensions, but a chain in higher 
dimensions cannot be reconstructed from the position of the vacancies. Moreover, if such a system is driven by a field, 
the transition rates not only depend on the vacancy distribution, but also on the direction of the connecting links. So 
one arrives at a model of the same complexity as the RD model. The model shows in d = 1 not the characteristics 
of the slow reptation behavior. The drift velocity for long chains approaches a constant rather than decaying with 
the inverse powers of the length. The reason is that the model has no obstacles which slow down the drift and the 
diffusion. 

The solution can also be used for a special RD model in d = 1, with the possibility of creation and annihilation 
of hernias. This shows the importance of the hernias in d = 1. Whereas it is believed that hernia creation and 
annihilation is of minor importance in higher dimension, it plays a decisive role in d = 1. Without the hernias as 
a move, the + taut links and the — taut links block each other and they are driven towards each other. A hernia 
annihilation followed by a creation of a pair of taut links in the opposite order, allows them to pass each other. So 
the hernias, which are very abundant in d = 1, effectively remove the obstacles, which are characteristic for the RD 
model. 
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